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ABSTRACT 

We report the discovery of the millisecond pulsar (MSP) PSR J1950+2414 (P = 4.3 ms) in a binary system 
with an eccentric (e = 0.08) 22-day orbit in Pulsar ALFA survey observations with the Arecibo telescope. Its 
companion star has a median mass of 0.3 Mg and is most likely a white dwarf. Fully recycled MSPs like this 
one are thought to be old neutron stars spun-up by mass transfer from a companion star. This process should 
circularize the orbit, as is observed for the vast majority of binary MSPs, which predominantly have orbital 
eccentricities e < 0.001. However, four recently discovered binary MSPs have orbits with 0.027 < e < 0.44; 

PSR J1950+2414 is the fifth such system to be discovered. The upper limits for its intrinsic spin period 
derivative and inferred surface magnetic field strength are comparable to those of the general MSP population. 

The large eccentricities are incompatible with the predictions of the standard recycling scenario: something 
unusual happened during their evolution. Proposed scenarios are a) initial evolution of the pulsar in a triple 
system which became dynamically unstable, b) origin in an exchange encounter in an environment with high 
stellar density, c) rotationally delayed accretion-induced collapse of a super-Chandrasekhar white dwarf, and 
d) dynamical interaction of the binary with a circumbinary disk. We compare the properties of all five known 
eccentric MSPs with the predictions of these formation channels. Future measurements of the masses and proper 
motion might allow us to firmly exclude some of the proposed formation scenarios. 

Subject headings: methods: data analysis, stars: neutron, pulsars: general, pulsars: individual (J1950+2414) 
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Millisecond pulsars (MSPs; Backer et al. 1982) are thought 
to be old neutron stars (NS) spun up by mass accretion and 
transfer of angular momentum from a companion star (Rad- 
hakrishnan & Srinivasan 1982; Alpar et al. 1982). The spin 
frequencies of these so-called “recycled” pulsars can range up 
to 716 Hz (Hessels et al. 2006; ATNF pulsar catalog, Manch¬ 
ester et al. 2005). The different evolutionary phases of these 
pulsars in binary systems and also the transitions between 
these phases have been recently observed in much more de¬ 
tail, though several major puzzles still remain. The binary 
system starts off as a low-mass X-ray binary in which a NS 
accretes matter from a companion star (Smarr & Blandford 
1976; Bildsten et al. 1997). The main emission from these sys¬ 
tem is X-rays from the hot accretion disk. These system can 
U'ansition into an accreting X-ray MSP in a binary (Wijnands 
& van der Klis 1998), in which matter is funnelled onto the 
neutron star’s surface and significant X-ray emission modu¬ 
lated by the NS spin is detected. After the accretion dies off, 
the NS can become “visible” as a radio MSP, powered by the 
rotation of the neutron star’s magnetic field (kauris & Savonije 
1999; Stairs 2004). In some cases, these systems are seen to 
switch on roughly year-long timescales between states as an 
LMXB and a radio MSP. For example, Archibald et al (2009) 
showed that the radio MSP PSR J1023+0038 has turned on af¬ 
ter a recent (~ 10 yrs) LMXB phase; more recently this system 
changed back into an LMXB (Archibald et al. 2014; Patruno 
et al. 2014; Stappers et al. 2014). Similarly, LMXB/radio MSP 
state transitions have been shown for PSRs J1824—2452 and 
J1227—4853 (Papitto et al. 2013; Bassa et al. 2014; Roy et al. 
2015). Together, these observations nicely demonstrate the 
basic recycling scenario; however, they have also raised new 
puzzles and shown that the details of the process are quite 
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complex. It is also possible that these three aforementioned 
transitional MSP systems, all “redback” MSPs (Roberts 2013), 
are not representative of the evolution of all types of MSPs 
and that these systems might not lead to recycled pulsars at all 
(Chen et al. 2013). As such, much remains to be understood in 
the formation of radio MSPs. 

The recycling pathway results in highly circular orbits of the 
binary system through tidal forces acting during the 10 s — 10 9 
year long accretion phase (Phinney & Kulkarni 1994). Until 
2008, this seemed true of all fully recycled MSPs. All binary 
pulsars with v > 50 Hz and v < 10 -14 Hzs -1 outside of 
globular clusters had orbital eccentricities between e = 10' ' 
and e = 10 3 (Manchester et al. 2005). These eccentricity 
limits do not apply to MSPs in globular clusters because their 
high stellar densities and resultant close stellar encounters can 
significantly increase the orbital eccentricity of a binary MSP 
after the end of the accretion phase (Rasio & Heggie 1995; 
Heggie & Rasio 1996) or even in some cases lead to exchange 
encounters. Indeed, several highly eccentric binary MSPs 
have been found in the Galactic globular cluster system (e.g., 
Ransom et al (2005); see also the online catalog of globular 
cluster pulsars at http : //www. naic . edu/ ~pf reire/ 
GCpsr. html). 

In 2008, the situation changed with the discovery of 
PSR J1903+0327. This fast-spinning MSP (P = 2.15ms) 
in an eccentric orbit (e = 0.44) with a ~ 1 Mq main-sequence 
star (Champion et al. 2008) cannot have formed through the 
“normal” binary evolution described above. Rather, it is be¬ 
lieved to have originated from a hierarchical triple that became 
dynamically unstable (Freire et al. 2011; Portegies Zwart et al. 
201 ; Pijloo et al. 2012). 

Soon afterwards, Bailes (2010) reported an “anomalous” 
orbital eccentricity of 0.027 for PSR J1618—3921 (Bailes; 
private communication), a MSP first reported in id wards 
& Bailes (2001). More recently, Deneva et al. (2013) and 
Barr et al. (2013) reported the discovery of two more un¬ 
usual binary MSP systems PSR J2234+06 (e = 0.13) and 
PSR J1946+3417 (e = 0.14). These three systems are fully 
recycled with spin periods between 3 and 12 ms, orbital peri¬ 
ods P or b from 22 to 32 days and median companion masses 
M 2 ~ 0.25 Mq; i.e., apart from the large orbital eccentricity 
all parameters are compatible with the canonical recycling 
formation channel leading to a MSP with a white dwarf (WD) 
companion as described above. 

These unusual orbital eccentricities require a non-standard 
formation channel. Like PSR J1903+0327, they could have 
formed in a triple system which later became unstable and 
ejected the outer (tertiary) companion. In fact, fully recycled 
MSPs can form in stable triple systems, as shown by the dis¬ 
covery of PSR J0337+1715, the first MSP in a stellar triple 
system with two white dwarf companions (Ransom et al. 2014; 
Tauris & van den Heuvel 2014). 

However, the chaotic disruption of a triple system would 
most likely not lead to the formation of MSPs with very similar 
orbital and spin characteristics as discussed above. These 
similarities suggest instead a more orderly mechanism with a 
more predictable outcome. 

A possibility for such a mechanism was proposed by reire 
& Tauris ( 1014). They suggest these systems formed through 
an accretion-induced collapse (AIC) of a super-Chandrasekhar 
mass oxygen-neon-magnesium white dwarf in a close binary. 
This star initially avoids AIC due to its rapid rotation. Only 
after the end of the accretion episode, and after the WD loses 
sufficient spin angular momentum, does it undergo AIC to 


directly produce an MSP in an eccentric orbit. A second 
possibility was suggested by Antoniadi; (2014): here the MSPs 
form through the usual channel, but the orbital eccentricity 
arises from the dynamical interaction with a circumbinary 
disk. This disk may form from donor material ejected during 
hydrogen-shell flash episodes. Antoniadi: (201 ) shows that 
even a short-lived disk can produce eccentricities as large as 
e = 0.15. 

Here we present the discovery and initial timing of the bi¬ 
nary MSP PSR J1950+2414 in Pulsar ALFA (PALFA) survey 
data obtained at 1.4 GHz with the Arecibo telescope. This 
is the fifth eccentric MSP in the Galactic field to be discov¬ 
ered; its orbital and spin parameters are similar to those of 
PSR J2234+06 and PSR 11946+3417. 

First, we will briefly describe the PALFA survey, the Ein¬ 
stein @Home project and its analysis of PALFA survey data, 
and the discovery of PSR J1950+2414. We then describe 
the timing observations, data reduction, and timing solution, 
followed by a discussion of PSR J1950+2414. We discuss 
our discovery in the context of different possible formation 
channels for MSPs in eccentric binaries. We conclude with 
an overview of future studies of this pulsar system and how 
they might allow us to exclude some of the possible formation 
channels. 

1.1. The PALFA Survey 

The PALFA Survey (Cordes et al. 2006) was proposed and 
is managed by the PALFA Consortium. It consists of about 
40 researchers (including students) at about ten institutions 
worldwide 21 . 

At Arecibo Observatory, the PALFA Consortium uses the 
Arecibo L-band Feed Array (ALFA 22 ). The output of the seven 
ALFA beams is fed into the Mock spectrometers 23 . The ob¬ 
serving band of 322.6 MHz is split into two overlapping bands 
with bandwidths of 172.0625 MHz each. The two sub-bands 
are centred at 1300.1680 MHz and 1450.1680 MHz, respec¬ 
tively. A total of 960 frequency channels is used, generated by 
polyphase filterbanks to enable the correction of radio pulse 
dispersion in the interstellar medium. Spectra are sampled 
every 65.4762 p,s. 

The PALFA survey setup (high observation frequency, large 
number of filterbank channels, fast-sampling spectrometers) is 
chosen to maximize the chances of discovering MSPs at large 
distances within the Galactic plane, where previous surveys 
have had little to no sensitivity to MSPs. At the same time, 
the detection of these objects has been difficult because of 
their high values of dispersion measure (DM), which induce 
high dispersive smearing per channel. The narrow channels 
used in the PALFA survey address this issue. Finding many 
MSPs is the highest priority in this survey because of their 
wide range of astrophysical applications: testing Einstein’s 
theory of general relativity and alternative theories of gravity 
(Antoniadis et al. 2013), measuring NS masses (Demorest 
et al. 2010), which can strongly constrain the equation of state 
of dense matter, finding suitable sources for pulsar timing 
arrays ( obbs et al. 2010), which will be used for detection of 
very low-frequency gravitational waves, improved estimates 
of the Galactic MSP population (Swiggum et al. 2014), and 

21 http : / /www2 . naic . edu/alfa/pulsar/ 

22 http : / /www . naic . edu/alfa/ 

23 Details of the Mock spectrometers may be found on the following 
NAIC web page: http : / /www . naic . edu/ -phil/hardware/pdev/ 
pdev . html 
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- in the case of PSR J1950+2414 and similar systems - a 
better understanding of stellar and NS formation and evolution. 
PAEFA discoveries are proof of the power of the survey to 
discover MSPs at high DMs (Champion et al. 2008; Knispel 
et al. 2010, 2011; Crawford et al. 2012; Deneva et al. 2012). 

Since its first observations in 2004, the PALFA Consortium 
has been surveying the part of the sky close to the Galactic 
plane (\b\ < 5°) that is visible to Arecibo Observatory, i.e., 
declinations 0° < S < 37°. The complete survey of this 
sky area will require about 330,000 separate Arecibo beams 
(equivalent to ~ 47,000 pointings of the seven-beam ALFA 
receiver). 

Data from the PALFA survey are analyzed by the Ein¬ 
stein® Home pipeline briefly described below and also by an 
independent pipeline operating on a supercomputer at McGill 
University using the PRESTO software package 24 (Lazarus et 
al., in prep.). 

1.2. Einstein@Home 

Einstein@Home 25 is a distributed volunteer computing 
project (Anderson et al. 2006). Members of the public do¬ 
nate otherwise unused compute cycles on their home and/or 
office PCs, and Android devices to the project to enable blind 
searches for unknown NS. 

Einslein@Home is one of the largest distributed computing 
projects. In the last ten years since its launch, more than 
390,000 volunteers have contributed to the project. On average, 
about 46,000 different volunteers donate computing time each 
week on roughly 105,000 different hosts 26 . The sustained 
computing power provided by these volunteers is currently of 
order 1.75 PFlop s _1 , which is ~5% the computing power of 
the world’s fastest supercomputer 

Einstein@Home analyzes data from the LIGO gravitational- 
wave detectors, the Fermi Gamma-ray Space Telescope, and 
large radio telescopes such as the Arecibo Observatory and the 
Parkes Radio Telescope. Observational data are stored and pre¬ 
pared for Einstein@Home processing on the Atlas computer 
cluster at the Albert Einstein Institute in Hannover, Germany 
(Aulbert & Fehrmann 2009). Volunteers’ computers download 
the data along with scientific software from dedicated servers 
in Hannover and Milwaukee, and run it automatically. The cen¬ 
tral coordination and management of the computing is handled 
by the Berkeley Open Infrastructure for Network Computing 
(BOINC; Anderson et al (2006)). 

The Einstein@Home PALFA analysis pipeline consists of 
three main steps. 1) Observational data are de-dispersed at 
~4,000 trial values to mitigate the radio pulse dispersion from 
the signal passing through the interstellar medium. Strong 
burst-like RFI is masked and periodic RFI is identified and 
replaced by random noise. Step 2) is done on the computers 
of the general public, attached to the project: each resulting 
de-dispersed time series is analyzed for periodic radio pulsar 
signals using Fourier methods. The Einstein@Home pipeline 
searches for radio pulsars in compact binary systems with 
orbital periods as short as 11 min. The required orbital demod¬ 
ulation to remove the Doppler effect from binary motion is 
done in the time domain. The demodulation is repeated for 
~ 7,000 different orbital configurations to cover a wide range 
of possible physical orbital parameters. For each of these, the 

24 https : //github . com/scottransom/presto 

25 http : / /einsteinathome . org 

26 http : // einsteinathome . org/server_status . html as of 

early April 2015 


Fourier analysis is repeated after the time-domain demodula¬ 
tion. The 100 statistically most significant candidates from 
each de-dispersed time series are stored and sent back to the 
project servers in Hannover. 3) The resulting ~ 400,000 candi¬ 
dates for each beam are sifted and remaining candidates folded 
using the raw data and graded by machine learning methods 
(Zhu et al. 201 ). Per-beam overview plots are also visually 
inspected to identify (of order a few) promising candidates in 
each beam. 

A full description of the Einstein@Home radio pulsar search 
pipeline employed to discover PSR J1950+2414 is beyond 
the scope of this publication and is available in Allen et al. 
(2013). The sifting techniques used to reduce the number of 
relevant candidates are described in Knispel (201 ); Knispel 
et al. (2013) 

To date, the project has discovered 51 neutron stars through 
their radio emission. As part of the PALFA collaboration, 
Einstein@Home has discovered a total of 27 radio pulsars 27 , 
including the fastest spinning disrupted recycled pulsar (i.e., a 
pulsar ejected from a binary system due to the companion star’s 
supernova explosion; Knispel et al. 2010), and a relativistic 
intermediate-mass binary pulsar (Knispel et al. 2011; Lazarus 
et al. 201 ). The Einstein@Home search for radio pulsars in 
archival data from the Parkes Radio Telescope has found 24 
pulsars missed by several previous re-analyses of this data set 
(Knispel et al. 2013). 

Einstein@Home has to date also discovered four gamma- 
ray pulsars in Fermi data ( : et al. 201 ). These ongo¬ 

ing searches use highly efficient data analysis methods ini¬ 
tially conceived for gravitational-wave data analysis (Pletsch 

& Allen 2009; Pletsch 2010, 2011; Pletsch et al. 2012). 

The detection of continuous gravitational waves from rotat¬ 
ing neutron stars in data from ground-based interferometric 
detectors is the main, long-term goal of Einstein@Home. The 
project has provided the most stringent upper limits for grav¬ 
itational waves from rapidly rotating neutron stars in blind 
searches to date (Abbott et al. 2009a, b; Aasi et al. 2013). 

1.3. Discovery of PSR J1950+2414 

PSR J1950+2414 was discovered by the Einstein@Home 
PALFA pipeline on 2011 October 4 with a statistical signif¬ 
icance S = 120.8 (negative decadic logarithm of the false- 
alarm probability in Gaussian noise) in a PALFA survey obser¬ 
vation from 2009 April 4. The pulsar was found at a spin period 
of 4.3 ms and a DM of 142 pc cm” 3 . The signal’s celestial na¬ 
ture was confirmed the same day with a second observation by 
the Arecibo telescope. 

2. OBSERVATIONS AND DATA ANALYSIS 
2.1. Observations 

Following its discovery, PSR 11950+2414 was observed 
during the PALFA survey observing sessions and with ded¬ 
icated observations by the 76-m Lovell telescope at Jodrell 
Bank. 

The observations with the Arecibo telescope employed the 
usual PALFA survey set up as described in Section 1.1, using 
only the central beam of ALFA. The durations were between 
4.5 mins and 15 mins. In total 24 PALFA survey observations 

27 All discoveries are available online at http:// 
einsteinathome . org/radiopulsar/html/rediscovery_ 
page/rediscoveries . html and http : //einsteinathome . 
org/radiopulsar/html/BRP4_discoveries/ 
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of PSR J1950+2414 were carried out from 2011 October 5 
until 2012 May 5. 

Between 2012 September 21 and 2012 December 31, nine 
dedicated timing observations of PSR J1950+2414 were car- + 
ried out as part of a PALFA millisecond pulsar timing cam- A 

paign. Observations used the L-wide receiver and the Puerto- c- 

Rican Ultimate Pulsar Processing Instrument (PUPPI) backend. 

A frequency band of 800 MHz in 2048 channels, centered on 
a frequency of 1380 MHz was observed in all sessions. The ~ 
usable bandwidth is limited by the L-wide receiver which de¬ 
livers an observation bandwidth of 700 MHz from 1.1 GHz to 
1.8GHz. Observation times were lOmins or 15mins with a 
sampling time of 40.96 /xs. 

The observations with the 76-m Lovell telescope used a 
cryogenically-cooled dual-polarization receiver at a central 
frequency of about 1520 MHz. A 512-MHz band was sam¬ 
pled at 8 -bits resolution and processed using a digital filter 
bank into 2048 0.25-MHz frequency channels. After radio 
frequency interference excision, approximately 384 MHz of 
usable bandwidth remained. The data were folded into 1024 
pulse phase bins and were de-dispersed, generating an aver¬ 
age profile for each 10-second sub-integration. Observations 
of mostly 30 mins duration were made approximately every 
10 days and the data set includes 61 TO As from the MJD period 
55851 to 56272. 

2.2. Data reduction 

Observational data were processed offline for radio fre¬ 
quency interference (RFI) mitigation, and to obtain times-of- 
arrival (TOAs). 

The separate sub-band data files obtained with the PALFA 
Mock spectrometers were merged into one file in PSRFITS 
format for each observation. PUPPI data were written in sepa¬ 
rate PSRFITS files, each covering adjacent observation time 
stretches. 

For RFI mitigation, the rfifind program from the 
PRESTO pulsar processing suite was used to generate an RFI 
mask for each observation. 

Each observation was folded at the appropriate dispersion 
measure and topocentric spin period with the prepfold pro¬ 
gram from PRESTO. Folded data files were used to create 
summed pulse profiles from all ALFA/Mock spectrometer ob¬ 
servations and all L-wide/PUPPI observations, respectively. 

Fig. 1 shows these ptdse profiles, averaged over 25 observa¬ 
tions with the Mock spectrometers and over ten observations 
with the PUPPI backend, respectively. The full width at half 
maximum duty cycle is 8 %. This corresponds to a width of the 
pulse of w 50 = 0.34 ms. The dispersive delay across a single 
frequency channel is 0.18 ms for the given DM, frequency res¬ 
olution and central frequency. An exponential scattering tail is 
apparent in the folded pulse profile. Centered on a rotational 
phase of 0.65 a small additional pulse component appears in 
both the ALFA/Mock and L-wide/PUPPI observations. 

We calibrate the profile using the radiometer equation fe.g., 
Lorimer & Kramer 2012) 


y/2A f T ohs /N 

to predict the observing system’s noise level o. We assume 
a system equivalent flux density 5 sys = 2.4 Jy for the L-wide 
receiver 28 . The combined observation bandwidth of L-wide 

28 Table 3 in http://www.naic.edu/~astro/guide/guide. 
pdf 



Rotational Phase 

Figure 1. The averaged radio flux density pulse profiles from the data col¬ 
lected with Arecibo Observatory using the Mock spectrometers and the ALFA 
receiver (dark-grey line) and with the PUPPI backend and the L-wide receiver 
(black line), respectively. The averaged radio flux density pulse profile from 
observations at the Lovell Telescope in Jodrell Bank is shown (light-grey line). 
The spin period of the pulsar is 4.3 ms. These profiles were used to obtain the 
TOAs. All pulse profiles have been normalized and aligned at their maximum. 
The period-averaged flux density at 1.4 GHz is .S 1400 — 119 zb 13 fiiy. 
with the PUPPI backend is A f = 700 MHz. We used a pulse 
profile summed from nine individual observations, six of which 
with 9.8 mins duration and three with 14.9 mins duration. The 
summed pulse profile was folded into N = 128 bins; then the 
expected off-pulse noise standard deviation used for calibration 
is 9.2 /xJy. Scaling the off-pulse noise level of the observations 
(between pulse phases 0.73 and 1.07 in Fig. 1), the resulting 
estimated period-averaged flux density of the pulsar at 1.4 GHz 
is 5*1400 = 119 + 13 /dy. 

The get_TOAs . py routine from PRESTO was used to ex¬ 
tract TOAs by employing the sums of pulse profiles as tem¬ 
plates. A single TOA was extracted for each ALFA/Mock 
spectrometer and L-wide/PUPPI observation, respectively. A 
total of 33 TOAs was extracted from the data obtained with the 
Arecibo telescope, 24 of which with the Mock spectrometer, 
the remaining 9 with the PUPPI backend. 

The PSRCHIVE 29 package by lotan et al (2004) was used 
for the data inspection, interference removal and arrival time 
determination in the Lovell Telescope data. A single TOA was 
generated for each of the observations using a high signal-to- 
noise template formed from the observations. 

The timing analysis used the TEMP02 software ( obbs et al. 
2006). Our timing solution uses the BT orbital model (Bland- 
ford & Teukolsky 1976; Edwards et al. 2006). The coherent 
solution fits 94 TOAs obtained between MJDs 55839 and 
56293, covering a baseline of 454 days. As visible in Fig. 2, 
the timing model accurately predicts TOAs. There are no clear 
trends in the TOA residuals as a function of time, nor as a func¬ 
tion of binary orbital phase. Constant time offsets between the 
three TOA data sets have been fitted with TEMP02. 

The weighted RMS of the timing residuals is 4.7 //s when 
including TOAs from both observatories with all receivers 
and backends. The most precise TOAs are obtained with the 
PUPPI backend and the L-wide receiver at Arecibo. For an ob¬ 
servation time of 15 mins, the TOAs have an average precision 
of 1.3 /xs. This makes PSR J1950+2414 an interesting pulsar 
to be used in Pulsar Timing Arrays ( obbs et al. 2010). 

3. TIMING RESULTS 

29 http : //psrchive . sourceforge . net 
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Figure 2. The TOA residuals of our timing model as a function of MJD 
(upper panel) and binary orbital phase (lower panel). Light-grey squares show 
TOA residuals from observations with the 76-m Lovell telescope, dark-gray 
circles residuals from observations with the Arecibo telescope using the ALFA 
receiver and the Mock spectrometers, and black triangles observations with 
the Arecibo telescope using the L-wide receiver and the PUPPI backend. The 
TOA residual distribution does not exhibit any clear trends. 
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Astrometric, spin, and binary parameters, including the rela¬ 
tivistic periastron advance, were determined from the timing 
analysis of the TOAs described in the previous section. Table 1 
shows the parameters of our timing solution obtained with 

TEMP02. 

The pulsar’s mass function is given by 

_ 4t T 2 xf _ Mf sin 3 (i) 

h T e P> b (M 1+ M 2 ) 2 ’ W 

where X\ = a\ sin (i) /c is the pulsar’s projected orbital 
semi-major axis in light-seconds, Tq = GMq/c 3 = 
4.925490947 fis is the solar mass in time units, P 0I b is the 
orbital period, and Mi and M 2 are the masses of the pulsar 
and the companion in units of solar masses, respectively. We 
find fi = 0.00626918(2) Mq, which indicates a low-mass 
companion with a minimum mass M 2 > 0.253 Mq, which 
is obtained for i = 90° and a pulsar mass M\ = 1.35 Mq. 
The median companion mass, assuming i = 60°, is M 2 , me A = 
0.297 Mq. Therefore, the binary companion is most likely a 
helium WD. 

Considering the system’s mass function, its high spin fre¬ 
quency v = 232.300 Hz and small spin period derivative, 
PSR J1950+2414 is very similar to the majority of comparable 
binary pulsar systems in the Galactic field. The most proba¬ 
ble nature for the system would therefore be a fully recycled 
millisecond pulsar with a WD companion. In this case, mass 
transfer from the WD progenitor to the MSP in the system’s 
past would have circularized the binary orbit, spun up the pul¬ 
sar, and damped its magnetic field, as described in Section 1. 

However, the orbital eccentricity e = 0.07981158(12) of 
PSR J1950+2414 is larger (by roughly two orders of magni¬ 
tude or more) than that of all millisecond pulsars with helium 
WD companions in the Galactic field with the exception of the 
three recently discovered similar systems (see Sections 1 and 

4). 

The theory of general relativity predicts an advance of the 
longitude of periastron, depending on the total mass of the 
system, its eccentricity and orbital period (Weisberg & Taylor 


Table 1 

Fitted and derived parameters for PSR J1950+2414. 


Parameter 

Value 0 

General Information 

MJD Range 

55839.0—56292.7 

Number of TOAs 

94 

Weighted RMS of Timing Residuals (fis) 

4.7 

Reduced-x 2 value 

1.2 

MJD of Period Determination 

55839 

Binary Model Used 

BT 

Fitted Parameters 

R.A., a (J2000) 

19:50:45.06390(10) 

Dec., S (J2000) 

+24:14:56.9638(11) 

MJD of Position Determination 

55839 

Spin Frequency, v (Hz) 

232.30014862462(14) 

Spin Frequency Derivative, z> (x 10 —15 Hzs -1 ) 

-1.020(6) 

Dispersion Measure, DM (pc cm -3 ) 

142.089(18)*' 

MJD of DM Determination 

55839 

Projected Semi-Major Axis, ai sin ( i ) (lt-s) 

14.2199738(11) 

Orbital Period, P or b (days) 

22.1913727(10) 

Epoch of periastron, To (MJD) 

55846.0226219(15) 

Longitude of periastron, cuo (deg) 

274.4155(3) 

Periastron advance, uj (degyr -1 ) 

0.0020(3) 

Orbital eccentricity, e 

0.07981158(12) 

Derived Parameters 

Spin Period, P (ms) 

4.304775549739(2) 

Spin Period Derivative, P (x 10 -20 ) 

1.8900(10) 

Intrinsic P, P int (xlO~ 20 ) 

<2.1166(10+ 

Galactic longitude, / (deg) 

61.10 

Galactic latitude, b (deg) 

-1.17 

Distance, d (NE2001, kpc) 

5 5 +0,8 

Mass Function, fi (Mq) 

0.00626918(2) 

Minimum companion mass M2 (Mq) 

0.253 

Total Mass, M (Mq) 

2.3(4) 

Characteristic Age, r c = P/ (2 P) (Gyr) 

>3.2 C 

Surface Magnetic Field Strength, Bg (x 10 s G) 

<3.1 c 

Spin-down Luminosity, E (x 10 34 ergs/s) 

< 1.0 C 


“ The numbers in parentheses are the \-a, TEMP02-reported uncertainties on 
the last digit(s). Uncertainties in the pulsar distance inferred from NE2001 are 
assuming a 20% uncertainty in DM to account for model uncertainties. 
b The solution here is based on single-frequency TOAs per Arecibo and 
Jodrell Bank observations, respectively. Fitting the DM is therefore degenerate 
with fitting the JUMPs between the data sets. The DM and its uncertainty 
were therefore fitted from TOAs obtained in two sub-bands for the Arecibo 
data combined with the Jodrell Bank TOAs, and kept fixed for the solution 
presented here. 

c The observed period derivative has been corrected for Galactic acceleration. 
We cannot correct for the Shklovskii effect, and therefore the value of -^mt is 
an upper limit, see Section 3. Its reported uncertainty here does not include 
the uncertainty from the Galactic model. r c , Bs, and E have been inferred 
from P ta t. 


198 ), given by 

" = 3 (++) -- ’ (3) 

where P or b, Mi, M 2 and 7+ are defined as above and e is the 
orbital eccentricity. 

Our observations of PSR J1950+2414 significantly detect 
the system’s periastron advance w = 0.0020(3) degyr -1 . This 
corresponds to an inferred total mass of M = Mi + M 2 = 
2.3(4) Mq, where the uncertainty is dominated by the uncer¬ 
tainty in lj. This measurement is suggestive of a NS signif¬ 
icantly more massive than 1.35 Mq, however, it is currently 
not precise enough to infer the individual masses of the pulsar 
and/or the companion star. We expect our ongoing observation 
campaign to significantly improve the measurement of ii) and 
therefore the estimate of the total mass. 

Ongoing timing observations might also help to constrain 
or measure the Shapiro delay in the binary system. Combined 
measurements of the Shapiro delay and the periastron advance 
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can be used to infer the separate masses of both components 

(Freire & Wex 2010). 

The observed spin period derivative P has been corrected for 
Galactic acceleration (Nice & Taylor 1995), to obtain a value 
of the spin period derivative of P con = 2.12 x 10~ 2 °. This 
takes into account a Galactic acceleration model that uses the 
radial velocity curve published in Reid et al. (2014). However, 
since the proper motion currently cannot be measured with 
useful accuracy, we cannot correct for the contribution Pshk 
to the period derivative from the Shklovskii effect (Shklovskii 
1970), and therefore we can only estimate an upper limit for 
the intrinsic Pi nt . 

Both the intrinsic Pj nt and the Shklovskii effect contribu¬ 
tion Pshk are always positive and fulfill Pi nt + Pshk = Peon- 
Therefore, assuming Pshk = 0, we obtain an upper limit to the 
intrinsic spin period derviative of l\ m < 2.12 x 1(T 20 . On 
the other hand, assuming P int = 0 we can constrain the total 
proper motion to < 19 mas/yr. 

From Pi n t = 2.12 x 10 20 s s -1 we infer a lower limit for 
the characteristic age of r c > 3.2 Gyr and upper limits for the 
the surface magnetic field strength Bs < 3.1 x 10 8 G and spin- 
down luminosity E < 1.0 x 10 34 ergs/s. These parameters are 
very similar to those of the general MSP population. 

3.1. Counterparts at other Wavelengths 

Let us now consider whether PSR J1950+2414 is observ¬ 
able at X-ray and 7 -ray wavelengths and whether the nature of 
its companion can be established by way of optical/IR obser¬ 
vations. 

The field around PSR J1950+2414 has not been previously 
targeted in a pointed observation by any former or current 
X-ray mission. A 581-s exposure from the ROSAT all-sky 
survey is too shallow to yield any useful constrains on the 
X-ray flux from the pulsar. Based on its E, PSR J1950+2414 
likely has a soft, thermal X-ray spectrum due to heating of its 
magnetic polar caps, with L x ~ 10~ 4 — 10 -3 P ~ 10 30-31 
erg s _1 , as typically seen in non-eclipsing MSPs (Zavlin 2006; 
Bogdanov et al. 2006). At a distance of 5.5 kpc and assuming 
an absorbing column of A+ «4x 10 21 cnT 2 (based on the 
empirical D.VI-A'h relation from He et al. 2013) a detection 
of the pulsar with Chandra or XMM-Newton would require 
impractically long exposures. 

There is no Fermi LAT Third Point Source Catalog (3FGL; 
The Fermi-LAT Collaboration 2015) object in the vicinity of 
PSR J1950+2414’s position. A visual inspection of the Fermi 
LAT data from the beginning of the mission up to 2015 January 
19 shows no excess emission at the pulsar position. This is ex¬ 
pected given that at the Galactic latitude of PSR J1950+2414 
(6 = — 1 . 1 °), the diffuse 7 -ray background is overwhelming. 
In addition, the commonly used 7 -ray “detectability” met¬ 
ric for pulsars, v E/D 2 is 3 x 10 15 erg 1 / 2 s -1 / 2 kpc -2 for 
PSR J1950+2414, substantially lower compared to MSPs de¬ 
tected by Fermi LAT (see, e.g.. Figure 15 in Abdo et al. 2013) 
so it is not expected to be a bright 7 -ray source. 

In principle, the nature of the companion of 
PSR J1950+2414 can be established by way of opti¬ 
cal/IR observations. If the secondary is a main sequence star 
- in this case PSR J1950+2414 would belong to the same 
class as of systems as PSR J1903+0327 - it is likely to be a 
M3.5 red dwarf based on the mass measurement from radio 
timing (~0.25 M 0 ). Data from the Digitized Sky Survey 
reveal no optical counterpart at the radio position of the MSP, 
likely because it is too distant and extincted. Indeed, using the 
relations by Predehl & Schmitt (1995) and He et al. (2013) 


combined with the pulsar’s dispersion measure, we estimate 
the visual extinction to the target to be Ay ss 2.3. Taking 
the absolute magnitude of a M3.5V star. My « 11.5 (Henry 
et al. 2002), and scaling with the dispersion measure derived 
distance («5.5 kpc) yields an expected apparent magnitude of 
V « 27.5 including extinction. 

If we consider the intrinsic colors of an M3.5 V star based on 
2MASS data (Lepine & Gaidos 201' ), and scale the extinction 
with wavelength (based on Cardelli et al. 1989), we obtain 
apparent magnitudes in the near-IR of J « 23.9, if ~ 23.1 
and, K « 22.6. Based on this, it is evident that near-IR imag¬ 
ing observations provide the better means of establishing the 
nature of the companion, which in turn helps constrain the evo¬ 
lutionary history of PSR J1950+2414. Existing near-IR data 
from the 2MASS all-sky survey (Skrutskie et al. 2006) show 
no counterpart at the pulsar position, although this provides 
no meaningful constraints on the PSR J1950+2414 compan¬ 
ion since the limiting magnitude is only ~16. Thus, deeper 
near-IR observations are required. 

If the companion star is a He WD, its absolute magnitude is 
15 > My >10 (Kaler 2006). The extinction given above and 
the DM-inferred pulsar distance yield an apparent magnitude 
of 31 > my > 26. Detection in the visual band in imaging 
observations with 8-meter class telescopes requires apparent 
magnitudes of my < 24 (Bates et al. 2015). Therefore, a de¬ 
tection or spectroscopic identification of an He WD companion 
in the optical band is currently impossible. 

4. DISCUSSION 

PSR J1950+2414 has an unusual combination of small pul¬ 
sar spin period and high orbital eccentricity. The small spin 
period and small spin period derivative (and inferred surface 
magnetic field) point to a long recycling episode in the sys¬ 
tem’s past, in which matter from the companion star accreted 
onto the pulsar, spun it up and reduced its magnetic field. Yet, 
the high orbital eccentricity appears at odds with this picture 
because long recycling episodes are believed to circularize the 
binary orbits to eccentricities in the range e = 10 -7 — 10" 3 
(Phinney 1992). 

Thus, these systems must form through different channels 
and/or experience some mechanism that increases the orbital 
eccentricity after the pulsar is recycled. We will discuss four 
proposed channels: 1) Formation in a hierarchical triple sys¬ 
tem which became unstable and ejected one of its members; 
2) Perturbations in regions with high stellar density (globular 
clusters); 3) The rotationally delayed accretion-induced col¬ 
lapse of a super-Chandrasekhar white dwarf; 4) Dynamical 
interaction of the binary with a circumbinary disk. 

Table 2 compares our discovery with the four known eccen¬ 
tric MSPs in the Galactic field. 

The similarity of the spin and orbital parameters 
of PSR J1950+2414 to those of PSR J2234+06 and 
PSR J1946+3417 (and also, to a lesser extent, to 
PSR J1618—3921) is striking, and is further evidence for 
the existence of a common formation mechanism that is 
clearly distinct from that of PSR J1903+0327 (which, as 
mentioned before, likely resulted from the chaotic disruption 
of a triple system, as indicated by its more massive main- 
sequence companion, much larger eccentricity and orbital pe¬ 
riod). PSR J1618—3921 differs from the aforementioned pul¬ 
sars by a longer spin period and less eccentric orbit, however 
the orbital period is similar to that of PSR J1950+2414, sup¬ 
porting the idea that something anomalous happens within this 
range of orbital periods. 
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Table 2 

Physical parameters of the five known eccentric MSPs. PSR J1903+0327 
(Champion et al. 2008) was most likely formed from a triple system. 
PSR J1618—3921 (data reproduced from Edwards & Bailes 2001, Bailes 
2010, and Bailes - private communication), two recently discovered MSPs 
(data from Freire & Tauris 2014; Barr et al. 2013; Deneva et al. 2013) and 
PSR J1950+2414 (from this publication) have very similar orbital and spin 
characteristics, suggesting a common formation channel. The height above 
the Galactic plane, z, was calculated from the DM-inferred distance and the 
Galactic latitude. 


Pulsar 

P 

Porb 

^2,med 

e 

1*1 

PSR J1903+0327 

2.1 ms 

95.2 days 

1.1M 0 

0.44 

0.1 kpc 

PSR J2234+06 

3.6 ms 

32 days 

0.23 Mq 

0.13 

0.4 kpc 

PSR J1946+3417 

3.2 ms 

27 days 

0.24 Mq 

0.14 

0.5 kpc 

PSR J1618—3921 

12.0 ms 

22.8 days 

0.20 M 0 

0.027 

0.6 kpc 

PSR J1950+2414 

4.3 ms 

22.2 days 

0.30 Mq 

0.08 

0.1 kpc 


All pulsars in Table 2 have small distances from the Galactic 
plane 1*1 < 0.6 kpc. If the systems were produced by 2) 
exchange interactions in globular clusters, their distances to 
the Galactic plane are expected to be large, similar to the values 
for the globular cluster population (Harris 1996). Therefore, 
it is unlikely that any of the pulsars originated in a globular 
cluster. Future observations of PSR J1950+2414 might be 
used to constrain its proper motion and further elucidate the 
issue. 

As described in Section 1, Freire & Tauris (201 ) recently 
have proposed the formation mechanism 3) for fully recycled 
MSPs in eccentric orbits. Since PSR J1950+2414 is such a 
system, let us now consider whether it could have also formed 
by the same mechanism. 

At the center of their model is an accretion-induced collapse 
of a super-Chandrasekhar mass white dwarf in a close binary. 
The WD only collapses after the end of the accretion episode, 
and after it has lost sufficient spin angular momentum, directly 
producing an MSP in an eccentric orbit. Based on their sim¬ 
ulations, Freire & Tauris (201 ) provide a list of predictions 
for MSP systems formed through the rotationally delayed AIC 
channel. At the current time, PSR J1950+2414 is consistent 
with these predictions as detailed in the following. 

From Fig. 3 in ire & Tauris (201), and assuming a 
range of pre-AIC WD masses between 1.37 and 1.48M 0 , 
isotropically directed kicks, and initial orbital periods between 
15 and 30 days, the kick velocity must be around 2kms -1 
to yield orbital eccentricities between 0.09 and 0.14. The 
orbital eccentricity of PSR J1950+2414 is just outside the 
lower end of that range (and that certainly is the case for 
PSR J1618—3921). This might suggest a slightly larger kick 
magnitude and a slightly wider range of eccentricities to be 
observed in these systems in the future, or a slightly smaller 
NS binding energy than that assumed in that simulation. 

The companion should be a He WD with a mass in the range 
0.24 — 0.31 M 0 , as predicted by the .P 01 -b-MwD relation from 
Fig. 5 in Tauris & Savonije (1999) for this orbital period of 22 
days. Although we cannot currently determine the nature of 
the companion of PSR J1950+2414, its median companion 
mass of 0.30 Mq is certainly compatible with this expectation. 

Freire & Tauris (201 ) predict a pulsar mass of typically 
1.22 — 1.31 Mq. Currently, we have not measured the total 
mass of the system to a sufficiently high precision to test this 
prediction. 

The system has low Galactic height, as expected from sys¬ 
tems that formed with very small kick velocity (< 10 km/s) in 
an AIC event. Indeed, the low Galactic heights of the MSPs 


in eccentric orbits discovered to date suggests small peculiar 
space velocities compared to the general MSP population. This 
assumption might be further tested with a future measurement 
of the proper motion of the system. 

The characteristics of PSR 11950+2414 are also consistent 
with dynamical interaction with a circumbinary disk (Anto- 
niadis 2014). The circumbinary disk model makes two predic¬ 
tions which could be tested for PSR J1950+2414 with further 
observations: 1) As in • re ire & Tauris (201- ), the compan¬ 
ions should also be He WDs with masses given by iauris & 
Savonije ( 999), which as we have seen agrees with our cur¬ 
rent mass constraint for the companion of PSR J1950+2414. 
2) The masses of the MSP and the peculiar space velocities 
should closely resemble those of circular binary MSPs, un¬ 
like the small velocities (< 10 km/s) predicted by the AIC 
model. With our current timing data, we cannot constrain the 
proper motion and can therefore neither confirm nor falsify 
this prediction. 

Whatever the formation scenario for systems like 
PSR J1950+2414 is, it has in this case produced an object 
with an inferred surface magnetic field comparable to the gen¬ 
eral MSP population. This does not falsify any of the models 
discussed above, but might be useful constraint on other possi¬ 
ble formation scenarios. 

The fact that the characteristic age (r c > 3.2 Gyr) is very 
similar to that of the general MSP population can be used to 
formulate a first rough estimate of the occurrence of these 
systems. It suggests that the relative frequency of occurrence 
of PSR J1950+2414-like MSPs in the population is compara¬ 
ble to the fraction of the currently known systems in the cur¬ 
rently known population. Assuming that PSRs J1618—3921, 
J1946+3417, J1950+2414, and J2234+06 all formed through 
the same mechanism and all have similar r c , the relative fre¬ 
quency of these systems is of order 4 of a total of 230 known 
MSPs 50 not associated with a globular cluster, roughly 2% of 
the population. 

5. CONCLUSIONS AND FUTURE WORK 

We have presented the discovery and initial timing of the 
fully-recycled (P = 4.3 ms) MSP PSR J1950+2414 in an 
eccentric (e = 0.08) orbit (P or b = 22 days) with a median 
companion mass of 0.3 Mq. PSR J1950+2414 is only the 
fifth system with a large eccentricity to be discovered. Its spin 
and orbital parameters are similar to those of three previously 
known systems with orbital periods in the range of ~ 20 — 
30 days, spin periods between 2 ms and 12 ms, large orbital 
eccentricities (0.03 < e < 0.15), and companions with masses 
0.2 Mq < m c < 0.3 Mq. 

This combination of parameters cannot be explained in the 
standard pulsar recycling scenario. The existence of now four 
known similar systems which probably have a helium WD 
companion suggests the existence of a common formation 
channel leading to these unusual MSP systems. The upper limit 
on the intrinsic spin period and the inferred surface magnetic 
field are comparable to those of the general MSP population. 

We compared the properties of PSR J1950+2414 with the 
predictions of four proposed formation channels. The initial 
evolution of the pulsar in a hierarchical triple and the origin in 
an exchange encounter in an high stellar density environment 
(e.g., globular cluster) are unlikely. Although we cannot con¬ 
clusively rule out or confirm formation through a rotationally- 

30 http : //astro .phys. wvu . edu/GalacticMSPs/ 

GalacticMSPs . txt 
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delayed AIC event or interaction with a circumbinary disk with 
the measurements presented here, future observations might 
allow us to test a variety of predictions made by the different 
possible formation models. 

The pulsar is currently being observed with Arecibo Obser¬ 
vatory, and this will allow us to better constrain the pulsar and 
companion masses, which is of highest importance to discrimi¬ 
nate between the AIC model and the circumbinary-disk model. 
The individual masses of the pulsar and the companion can 
be obtained from combining observations of the relativistic 
periastron advance with observations of the Shapiro delay. Us¬ 
ing the orthometric parametrization of the Shapiro delay from 
Freire & Wex (2010) should provide the most precise estimates 
of the individual masses. The uncertainties of the periastron 
advance and the Shapiro delay in this parametrization are less 
degenerate than in the standard parametrization. Therefore, 
even if the measurement of the Shapiro delay has larger errors, 
a precise measurement of the relativistic periastron advance 
will allow for well-defined pulsar mass determination, see, e.g.. 
Lynch et al (2012). The ongoing observations will improve 
the relativistic periastron advance measurement and might also 
lead to a detection of the Shapiro delay. 

Radio observations with longer time baselines will pro¬ 
vide upper limits on or a measurement of the proper mo¬ 
tion of the pulsar, and therefore of its peculiar space velocity. 
The AIC model predicts small values (< 10 km/s), while the 
circumbinary-disk model predicts larger values. Measuring the 
proper motion could therefore permit discrimination between 
the two models. 

Further information about the nature of the companion might 
be obtained from observations at optical or infrared wave¬ 
lengths, as detailed in Section 3.1. Near-IR imaging obser¬ 
vations provide the best means of establishing the nature of 
the companion, which in turn helps constrain the evolutionary 
history of PSR J1950+2414. 

Modem pulsar surveys like PALFA are probing the Galactic 
disk to unprecedented depths and are likely to find other such 
systems. Having a larger sample size will allow discriminating 
between the different proposed models with greater confidence 
than is possible now. 
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